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I. INTRODUCTION 



A commutated network is an electrical network, either active or 
passive, consisting of resistors and reactive elements which are 
switched or commutated in a periodic manner. In most practical cases 
the reactive elements are capacitors. There are several configurations 
of a commutated network used in the field of electrical engineering. 

The configuration used here is that of commutated capacitors used in 
the feed-back path around an operational amplifier. Thus by applying 
commutation action to the capacitors in a lag network, the low-pass 
frequency response of that network can be transformed into that of a 
highly selective bandpass filter. The center frequency of the bandpass 
characteristic will be the same as the commutation frequency. 
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II. COMMUTATED NETWORK CONFIGURATION AND 
COMPUTER SIMULATIONS OF THE NETWORK 



The fundamental element of commutated networks is the capacitor 
shown in fig. 2.1. C 
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Figure 2.1 Commutated capacitor 
The commutating function, P, is a unit square wave of the same 
frequency as the input signal, e ^ . The following equations can be 



written for fig. 2.1: 






i dt 
c 



i = P* i . 
c 1 

e = P* e 



- p - i / p - 1 1 
■ ih « 



dt 



(2-1) 

( 2 - 2 ) 

(2-3) 

(2-4) 

(2-5) 



since P-P = 1 both when P is in the low and in the high position. 

Figure 2.2 shows how the capacitor is placed in the coupled RC 
commutated network. I|C 
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Figure 2.2 Coupled commutated network 
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Another variation of the commutated network is the uncoupled RC network 



shown in fig. 




-o 




Figure 2.3 Uncoupled commutated network 
Reference 1 proves that both variations of the commutated network 
can be represented by the circuit shown in fig. 2.4. 




Figure 2.4 Equivalent circuit of commutated networks 
The number of capacitors in the commutated network can be varied. 
For the following, four capacitors are used in the coupled version of 
the network. Figure 2.5 shows the circuit used in the analysis. 




Figure 2.5 Four-capacitor commutated network 
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The equivalent network for the four-capacitor coupled commutated network 



is shown in fig. 2.6. 




A computer program was written - in IBM's Digital Simulation 



Language - to simulate the network of fig. 2.5. Figure 2.7 shows the 



program used. 



F = 1. 

X = SIN(F*2.*3. 1415*TIME) 

ZONE = IMPULS (.0,1.00000) 

ZTWO = IMPULS (.1250,1.00000) 

Z THREE = IMPULS (.2500,1.00000) 
ZFOUR = IMPULS (.3750,1.00000) 
FONE = PULSE (ZONE, .499) 

FTWO = PULSE(ZTWO, .499) 

FTHREE = PULSE (ZTHREE, .499) 
FFOUR = PULSE(ZFOUR, .499) 

PONE = 2 . *FONE -1. 

PTWO = 2 . *FTWO -1. 

PTHREE = 2 . *FTHREE -1. 

PFOUR = 2 . *FFOUR -1. 

XONE = -X - . 95*Y 
WONE = PONE*XONE 
WTWO = PTWO*XONE 
WTHREE = PTHREE*XONE 
WFOUR = PFOUR*XONE 
OONE = INTGRL ( . 0 , WONE ) 

OTWO = INTGRL( . 0 ,WTWO) 

OTHREE = INTGRL( . 0 ,WHTREE) 

OFOUR = INTGRL(.0, WFOUR) 
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YONE = PONE*(-OONE) 

YTWO = PTWO*(-OTWO) 

YTHREE = PTHREE * ( - OTHRE E ) 

YFOUR = PFOUR* (-OFOUR) 

Y = - ( YONE + YTWO + YTHREE + YFOUR) 

Figure 2.7 Computer program for simulation 
The optimal value of the feed-back gain, K, was used. Reference 2 
proves that the optimal value is 0.95 ( 1/ ( 1H-1/7 +1/9 +1/15 +1/17 + . ..). 

The modulating square wave, p, was generated in the following way: 
an impluse, Z^, was generated at the beginning of every period. This 
impulse triggered a square-wave generator and the result, F_^ , is shown 
in fig, 2.8. 
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Figure 2.8 Square wave, F^ 

The square wave, F^ , was multiplied by a factor of two and the DC-level 

was shifted minus one unit. The result is P. as shown in fig. 2.9. 
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Figure 2.9 Square wave, P 

The different square waves are delayed with respect to one another. P^ 
is delayed T/8 with respect to P^. P^ is delayed T/8 with respect to P^ , 
and P^ is delayed T/8 with respect to P^ . t^ was chosen equal zero. 
Figure 2.10 shows the square waves as used in the program. 
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The following equations can be written for the commutating waveforms. 

is a unit square wave of period T, starting at t=0 thus f = 1/T Hz 

and cd = 2tt/T. 
n 







Figure 2.10 Commutating square waves 
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(2-6) 




(2-7) 



( 2 - 8 ) 



The objective of the RC-commutated network is to construct a 
signal of same frequency but of opposite phase as the input signal. 

If this is the case the sum of the input signal and the output signal 
has zero component of the input frequency. 

Figure 2.11 shows how this can be implemented. 



Here OUT is the sum of X and Y. As a result of the square waves being 
applied to the capacitors, harmonics of the frequency of the input 
signal are present in the signal OUT. If N is the number of commutating 
capacitors in the network, ref, 2 shows that the following frequencies 
are present in the signals OUT and Y: (2N-l)co n , (2N+l)cc> n , (4N-l)co n , 

(4N+1)cjo , etc. 



In order to analyze the discontinous signal Y, a computer program 
was written. The object of the program was to calculate the components 
a and b in a series of the form 



X 




OT 



2.11 



Figure 2.11 Filter configuration 
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(2-9) 
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Figure 2.12 shows the program written in FORTRAN IV. Table I lists 
the component values a^ and of the output signal Y. 

It is seen from the table that the following components are dominant: 
1st, 7th, 9th, 15th, 17th, 23rd, 25th, 31st, 33rd, 39th, 41st, 47th, 49th, 
etc. This is in agreement with theory given in ref. 2. 

An investigation was made to find out if the value of t^, the time 
where first goes to the value +1.0, had any effect on the output 
signal Y. The following phase delays were tried for the commutating 
wave - forms : tt/5, 2tt/5 , 3tt/5, and 4rr/5. The results showed that phase- 

delaying the commutating square waves does not change the output signal. 
The series of equation (2-9) still has the same coefficients. 
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FORTRAN 

SUBROUTINE TORIT(FNT,N,M,A,B) 

DIMENSION A(100) ,B(100) ,FNT(10500) ,C(100) ,PHI(100) 

WRITE (6,11) 

11 FORMAT (7X, ' J ' , 8X, 'A(J) ' , 13X,'B(J)', 13X,'C(J)' , 13X, 

60 AN=N 

COEF=2 . 0/ (2 . 0*AN+1 . 0) 

CONST=3 . 141593*COEF 
S l=SIN(CONST) 

Cl=COS (CONST) 

P = 1.0 

S=0..0 

J=1 

FNTZ=FNT(1) 

70 U2=0 . 0 
U1=0.0 
1=2 *N+ 1 

75 UO=FNT(I)+2 . 0*P*U1-U2 
U2=U 1 
U1=U0 
1= 1-1 

IF(I-l) 80,80,75 
80 A(J)=C0EF*(FNTZ+P*U1-U2) 

B(J)i=COEF*S*Ul 
K = J-l 

C (J) = SQRT((A(J))*(A(J) ) + (B(J))* (B(J))) 

PHI (J) = -ATAN(B(J)/A(J))*57 .2957795 
WRITE (6, 12) K,A(J),B(J),C(J),PHI(J) 

12 FORMAT (3X, 15, 1P4E 17 .7) 

IF (J- (M+l) ) 90,100,100 

90 Q=C1*P-S 1*S 
S=C1*S+S1*P 
P = Q 
J=J+1 
GO TO 70 

100 A(1)=A(1)*0.5 
RETURN 
END 



Figure 2.12 Computer program for Fourier analysis 



'PHI(J)') 
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N A B 



0 


-4.0915914E-02 


0.0 


1 


5 . 7432428E-02 


-9 . 7068429E-01 


2 


-6 .439543 8E-04 


-2.9122231E-03 


3 


6 . 73 18883E-04 


-3.3683525E-03 


4 


-8.5144676E-04 


-6. 864587 8E-04 


5 


-1.0724834E-03 


-1.2038115E-03 


6 


-1.3770524E-03 


4.9475231E-05 


7 


1 . 2983000E-01 


-5 . 76 17463E-02 


8 


-1.4599903E-04 


-4.4301525E-04 


9 


9.9531114E-02 


-5.0481014E-02 


10 


1. 1894142E-03 


-1.2319605E'-03 


11 


7.7437004E-04 


-5.9847650E-04 


12 


2. 1180733E-04 


-9. 1678812E-04 


13 


2.8510741E-04 


- 1 . 2260592E-04 


14 


-5.2539515E-04 


-1.4881822E-03 


15 


4.4528384E-02 


4 . 89047 10E-02 


16 " 


3 .3075595E-04 


-3.5192422E-04 


17 . 


4.2879425E-02 


4. 1087449E-02 


18 


1.3891549E-03 


5 . 6604668E-04 


19 


7 . 9694903E-04 


3 . 0924985E-04 


20 


8.3581335E-04 


-1.5659863E-04 


21 


1 .03853 11E-04 


-2.5742664E-04 


22 


1 . 0932467E-04 


-9 . 732847 IE-04 


23 


-1. 9760877E-02 


3 . 81517 86E-02 


24 


4 . 1980087E-04 


1 . 3531574E-04 


25 


- 1 . 6020168E-02 


3 . 723 1848E-02 


26 


-5 . 6435747E-05 


1.4239037E-03 


27 


5.6750767E-05 


8.6518889E-04 


28 


4.2703399E-04 


7. 1566226E-04 


29 


6 . 2502408E-04 


8. 1200222E-04 


30' 


1.3039294E-03 


6.7783 101E-04 


31 


-3 . 1548135E-02 


-3 . 6806425E-03 


32 .. 


3. 1624935E-05 


4 . 2581698E-04 


33 


-3.0806191E-02 


-1.3784338E-03 


34 


-1.3438673E-03 


4.3319259E-04 


35 


-7 . 84153 12E-04 


3 . 8454239E-04 


36 


-5 . 2506500E-04 


6.4045424E-04 


37 


-5 . 08794 10E-04 


8 . 5357740E-04 


38 


-1.7 829986E-04 


1.4619438E-03 


39 


-6. 1985478E-03 


-2.4380792E-02 


40 


-3 . 8792 193E-04 


1 . 8413758E-04 


41 


-7 . 8365095E-03 


-2.3663379E-02 


42 


-8. 8027469E-04 


- 1 . 1022575E-03 


43 


-6 . 67 1585 IE-04 


- 5 . 937633 IE -04 


44 


-7.9095364E-04 


-2 . 7252291E-04 


45 


-9 . 7992225E-04 


-1. 553283 IE-04 


46 


-1.4365583E-03 


3 .47 10323E-04 


47 


1 . 697 8830E-02 


- 1 . 1989255E-02 


48 


-3 . 2683252E-04 


-2.9340992E-04 


49 


1 . 62 12072E-02 


-1.3255987E-02 


50 


7 . 035.7 858E-04 


-1.2116958E-03 



Table I Fourier coefficients of output signal Y 
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Ill, FREQUENCY RESPONSE OF FILTER 



An investigation of the phase difference between output and input 
signals was done for frequencies within one decade on either side of 
the notch frequency, o^. 

The frequency response was done for a notch frequency of one hertz 
and with a time constant of one second. 

The equivalent network of fig. 2.6 was used for the frequency re- 
sponse study where the frequency of the input signal, x, was varied 
between 0.1Hz and 10.0 Hz. 

A Fourier analysis was done of theoutput signal y. The program of 
figure' 2.12 was used for the analysis. The phase difference between 
the input signal and the first harmonic of the output signal was calcu- 
lated. The calculations were done for times when the phase of the 
input signal was zero. The phase, of the first harmonic of the 

output signal is 

- tan (B^/A^) (3-1) 

Figure 3.1 shows the plot of phase versus frequency as obtained. 

The phase starts at zero degrees for low frequencies and goes to -90 
degrees for frequencies slightly less than jumps to +90 degrees for 

frequencies slightly above a) , and goes back to zero degrees for high 
frequencies . 

The values obtained in the computer simulation are collected in 
Table II. 
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Figure 3.1 Frequency response of filter - phase 
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Frequency 


Phase 


0.1 


- 1.51 


0.2 


- 7.40 


0.25 


-11.90 


0.3 


-18.59 


0.4 


-20.12 


0.5 


-24.00 


0.6 


-30.07 


0.8 


-51.54 


0.9 


-71.23 


1.25 


43.80 


2.0 


15.14 


2.5 


14.33 


4.0 


4.38 


10.0 


2.09 


20.0 


0.28 



Table II Results of frequency response simulation - phase 
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I V . A PPLI CATIONS OF FILTER IN S E RVOSYSTE MS 



Computer simulations were done with the filter used in conjunction 
with servosys terns . The objective of the filter was to suppress the 
component of frequency go^ in the output signal C. 

Example 1 

A servosystem with the following transfer function was used in the 
first example: 



c(s) 

xl(s) 



KMO/co 2 

n 

2 

s (0.00025 s+1) (^— j s+l) 

go n 



(4-1) 



n 

Figure 4.1 shows the system that has the open-loop transfer function 
of equation (4-1) . 




Figure 4.1 Servosystem of example 1 

The following values were used: ^ = 0.01, go^ = 2iTf^ 2500 = 15708 rad/sec. 

Figure 4.2 shows how the filter was inserted in the servosystem. Figure 
4.3 shows the Bode-plot for the system. 
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4.2 




Figure 4.2 Compensated filter network 
The signal VN is the sum of the signals OUT and VM. Since the funda- 
mental component (frequency = cjd^) of the signal OUT is 180 degrees out 
of phase with the signal VM, the sum signal VN will have no component 
of frequency cjd^. The computer simulation of the system of fig. 4.1 was 
carried for different values of the gain KMO. Figures 4.4 - 4.6 show 
the time responses of the system to a unit step input for different 
values of KMO. It is seen that the output of the system when the 
filter is inserted is much smoother than the output of the system 
without the filter inserted, and it does not contain any component of 
frequency cjd^. 

It is also seen that as the gain KMO is increased, the oscillations 
of the filtered output also increase in magnitude, but their amplitudes 
in the filtered output are still much smaller than their amplitudes in 
the unfiltered output signal. While the unfiltered output signal be- 
comes unstable, the filtered output is still stable. An analysis of 
the output signal C shows that the dominant frequency of the signal 
varies as a function of the gain KMO. The dominant frequency is found 
to be the frequency where the magnitude curve of the Bode-plot crosses 
the zero-decibel axis. 



25 



loo. tooo. toooO- 




26 



Figure 4.3 Bode plot of system of figure 



Example 2-4 

In examples 2-4 the same servosystem was used. Only the value of the 
gain KMO was changed. Figures 4.7 - 4.9 show the Bode-plots for the 
servosystem with and without the filter inserted. In example 2 the 
notch frequency is set at the phase cross-over (frequency = 15708 rad/sec). 
In example 3 the notch frequency of the filter is set at the magnitude 
cross-over of the Bode-plot. Finally in example 4 the notch frequency 
is set one decade below magnitude cross-over (frequency = 314 rad/sec). 

The computer simulations for the time responses for examples 2-4 
agreed with the expected results. 

In example 2 the unfiltered output is stable with a phase margin of 
about 13 degrees. The output when the filtered is in the network has 
less phase margin but is still stable. 

In example 3 the unfiltered output is stable with a phase margin 
of about 18 degrees. As seen from the Bode-plot for the system with 
the filter in the circuit the output is now unstable. 

In example 4 the unfiltered output is identical to that of example 
3. The notch frequency, however, is now one decade below the frequency 
a)^, and the filter's frequency responses are now back to Odb and 0 
degrees for magnitude and phase plots respectively at the gain cross- 
over. The filter has thus no influence on the output signal. 
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Figure 4.4 Time response to unit step input.KMO=8.636xl0 11 
Scale: X=0.0004 sec/in, Y=0.2 units/in. 
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12 

Figure 4.5 Time response to unit step input .KMO^ . 0x10 
Scale: X=0.0004 sec/in, Y=0.2 units/in. 
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Figure 4.6 Time response to unit step input ,KM0=2 . 5x10 
Scale: X=*0.0004 sec/in, Y=*0.2 units/in. 
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Figure 4.7 Bode-plot for example 
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Figure 4.8 Bode-plot for example 



loooo 
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Figure 4.9 Bode-plot for example 



V. APPROXIMATION OF FILTER TRANSFER FUNCTION 



As seen from fig. 3.1 the phase part of the frequency response is 
discontinous at the notch- frequency , 

A transfer function of the form 



is seen to have the same phase characteristic. By varying the value 
of the damping factor the curvature of the phase curve can be changed. 
Figure 5.1 shows the locations of the poles and zeros of the transfer 
function T(s) in the s-plane. The poles are shown for different values 
of the damping factor. 




s + a) s + a) 



(5-1) 



n n 




-§= 0.1 




§=0.3 



n 






Figure 5.1 Poles and zeros of T(s) 
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Figure 

filter 



5.2 shows the Bode-plot of T(s) and of the commutating network 
as obtained in chapter III. 
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Figure 5.2 Bode-plot of T(s) - phase 



VI. CONCLUSION 



The commutated network filter has a highly selective bandpass 
filter characteristic where the center frequency of the bandpass 
frequency-response curve can be easily varied by changing the frequency 
of the commutating functions. 

The filter can be approximated by a linear second-order transfer 
function for design purposes. The transfer function has two zeros on 
or very close to the imaginary axis and two poles on a circle with 
radius The location of the poles can be changed depending on the 

value of the damping factor. 
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APPENDIX A 



TIME RESPONSES OF SERVOS YSTEM 

A compensated servosystem with open-loop transfer function 

KMO 



c(s) _ 
xi(s) 



(as+1) (s^+2§ cjd s4a) 2) 
3 n a 



(A- 1 ) 



was used in example 1, tihapter IV to get time responses to a unit step 

input for the system with and without the filter inserted. 

Additional simulations were carried out for different combinations 

of the natural frequency cd^ and the damping factor Figure A.l shows 

which combinations of (jd and § were used. 

n 3 



.015 



.025 





f=2600 

f=2550 

f=2500 

f=2450 

f=2400 



f=2250 



<s 



Figure A.l Values used for time responses 
Table III lists the values used for KMO, (jd and 
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Figure 


CD /2tt 
n 


5 


KMO/10 11 


X-sca le 
sec/ in 


Y-scale 
units/ in 


A. 2 


2450 


0.010 


8.813 


0.0004 


0.2 


A. 3 


2550 


0.01 


8.463 


0.0004 


0.2 


A. 4 


2500 


0.005 


8.636 


0.0004 


0.2 


A. 5 


2500 


0.015 


8.636 


0 . 0004 


0.2 


A. 6 


2400 


0.01 


8.986 


0.0004 


0.2 


A. 7 


2600 


0.01 


8.290 


0.0004 


0.2 


A. 8 


2500 


0.0025 


8.636 


0.0004 


0.2 


A. 9 


2500 


0.020 


25.000 


0.0004 


0.5 


A. 10 


2500 


0.020 


8.636 


0.0004 


0.2 


A. 11 


2250 


0.010 


9.499 


0.0004 


0.2 


A. 12 


2250 


0.0025 


9.499 


0.0004 


0.2 


A. 13 


2250 


0.020 


9.499 


0.0004 


0.2 


A. 14 


2750 


0.01 


7.772 


0.0004 


0.2 


A. 15 


2750 


0.0025 


7.772 


0.0004 


0.2 


A. 16 


2750 


0.02 


7.772 


0.0004 


0.2 



Table III Values used for computer-simulation 
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Figure A. 2 Time response to unit step input 
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Figure A. 3 Time response to unit step input 
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Figure A. 4 Time response to unit step input 
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Figure A. 5 Time response to unit step input 
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Figure A. 6 Time response to unit step input 
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Figure A. 7 Time response to unit step input 
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Figure A. 8 Time response to unit step input 
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Figure A. 9 Time response to unit step input 
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Figure A. 10 Time response to unit step input 
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Figure A. 11 Time response to unit step input 
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Figure A. 12 Time response to unit step input 
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Figure A. 13 Time response to unit step input 



51 




Figure A. 14 Time response to unit step input 
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Figure A. 15 Time response to unit step input 
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Figure A. 16 Time response to unit step input 
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